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Abstract: A chaotic system under periodic forcing can develop a periodically 
visited strange attractor. We discuss simple models in which the phenomenon, 
quite easy to see in numerical simulations, can be completely studied analytically. 
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In chaotic systems synchronization phenomena, i.e. resonances, can 
occur: the simplest being motions asymptotically developing over a low 
dimension attractor which is visited in synchrony with a weak time- 
periodic perturbing force. This is analogous to the synchronization 
("phase locking") occurring when an integrable system, perturbed by a 
periodic forcing, develops attracting periodic orbits with period being 
a simple fraction (1, 3/2, 2,. . . ) of the forcing period. Formation of an 
attractor and existence of dissipation, i.e. of phase volume contraction, 
are concomitant. A chaotic system under a weak volume preserving 
\0 ' periodic perturbation can have, instead, a very singular chaotic evolu- 

tion, [1, 2], with dense orbits. A simple illustration of the resonance 
phenomenon in a chaotic system is exhibited here and exactly solved 
by deriving the shape of the attractor: which turns out to be a continu- 
ous surface with weak smoothness properties (being a strange attractor 
with an exponent of Holder-continuity with a very small lower bound). 
The peculiarity of the class of systems treated is the possibility of de- 
termining analytically an attractor whose existence can be rather easily 
^ \ seen in simulations. 

M. 

1 Resonances and chaos 

A chaotic system weakly interacting with a very regular system may generate a 
fully chaotic dynamics: as a basic example, consider the suspension flow of Arnold's 

cat map A2 = ( ^ j ^ with a flat floor function (i.e. the flow on T 3 = [0, 2ir} 3 

where every point x £ T 2 of the base moves vertically with constant velocity until 
it reaches the flat floor, at which point it reappears on the base, at the location 
A2X) coupled with a clock, which moves at the same speed as the suspension flow. 
If observed at the times when the suspension flow hits the floor, the unperturbed 
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dynamics takes the form of a linear discrete dynamical system A3 : T 3 — > T 3 , with 
if t 2 1 °\ 

A 3 = 1 1 . Consider now a small perturbation of A3: it is remarkable that 

Vo 1/ 

in the case of volume preserving perturbations, some pathological phenomena can 
take place [1, 2]. Namely, there is an open set U of C 1 conservative perturbations 
of A3 such that the following happens: the perturbed motion g G U is ergodic with 
respect to the volume measure in T 3 and there exists a foliation of T 3 into C 2 curves 
l x that are invariant under g (i.e. gt x = i gx )', moreover, the central Lyapunov 
exponent along the leaves is positive on the set $7 of full three-dimensional measure 
where it is defined; the set f2 intersects every leaf at exactly k points, for some 
k > [1, 2]. Such results provide an interesting example of an absolutely singular 
foliation, in the sense of Katok [3]. 

The situation is destined to change substantially in the presence of friction: if 
we couple via a dissipative perturbation the suspension flow of A2 with a clock, 
the phase between the two (that is arbitrary in the unperturbed case) might lock, 
making the motion fall onto a lower dimensional attractor, which would naturally 
be close to the union of k elements of a "horizontal" foliation of T 3 . We can 
say, when this happens, that there is a resonance between the periodic forcing 
and the chaotic evolution on T 3 : and unlike the cases quite typical in generic 
perturbations of periodically forced integrable systems where the resonance may 
lead to the appearance of stable periodic orbits, in the presence of hyperbolicity 
and of small dissipation the forcing gives rise to a "periodic strange attractor", 
i.e. to an attractor for the time evolution observed at times integer multiples of 
the forcing period, on which the dynamics is conjugate to a uniformly hyperbolic 
system. 

In this paper we provide examples of the above phenomenon: remark that 
in the non-dissipative case the perturbations discussed in [1, 2] can be in a set 
arbitrarily close to the unperturbed system and open in the C 1 topology (not 
containing the unperturbed system itself). Here we consider a class of analytic 
perturbations depending on a parameter e and show that for any small £ / 
the resonance occurs. The phenomenon is reminiscent of certain synchronization 
phenomena in chaotic systems [4] and the simple case considered below is, as far 
as we know, one of the first explicit examples of such phenomena that can be fully 
worked out analytically. 
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2 A simple example 

Denoting by x = (xi,X2) a point of T 2 , consider the evolution equation: 

x = 5(z) (Sx — x) + e f(x, z, w), 

w = 1 + e g(x, z, w), (2-1) 
z = l, 

with (x, ui,z)eT 2 xTxT and e > 0. Here 5 is the 2-7r-periodic delta function, 
with the convention that Jq dt' 6(t') = 1 for all t > 0. Set / = to simplify as much 
as possible the problem. The map S could be any linear hyperbolic automorphism 
of T 2 , in particular S = Ai- 

Then we can prove the following result. 

Theorem: Let S* be the flow on T 4 associated with Eq.(2.1) with f = and g 
analytic. Let also 

j Q (x,w) = f J** dtg(Sx,t,w+t), 7i0c,«0 = f dt -^-g(Sx,t,w+t), (2.2) 

and suppose existence of wo such that 7 (x,wo) = and ji(x,wo) = T < 0, 
independent of x,w. Then there are constants Eo,C > such that for < e < £q 
there exists a Holder- continuous function W : T 2 — >■ T of exponent (3 > Ce and 

• the surface A d = {(x, W(x), 0) : x G T 2 } is invariant under the Poincare 
map Sg W ; 

• the dynamics of Sg n on A is conjugated to that of S on T 2 , i.e. 

S^ix, W(x),0) = (Sx, W{Sx),0) . (2.3) 

Remarks. 

(1) It should be clear from the proof below that, for e > 0, A is an attractor for 
Sg with negative central Lyapunov exponent equal to |^(l + 0(e)), see concluding 
remarks in Section 7. Furthermore the function W is holomorphic in e = pe ta 
inside the domain which in polar coordinates is p(a) = £o cos a. 

(2) A simple explicit example of a function g satisfying the assumptions of the 
theorem with wq = ir and T = — 1 is: 



g(x, z, w) = sin (w — z) + sin (x\ + w + z) 



(2.4) 
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(3) It would be interesting to estimate the positions of the poles of the power 
spectrum of the correlations of pairs of smooth observables in the sense of the 
chaotic resonances theory in [5]. The spectrum is meromorphic in a strip with 
poles off the real axis for a dense set of observables: the unperturbed system has 
poles at — ilogA + + 2irk, k G Z, if A + is the largest eigenvalue of the matrix S; 
we expect the position of the poles closest to the origin to be modified by order e. 
Notice that in the unperturbed system the spectrum is an entire function for all 
analytic observables and the poles really appear with non-zero residue for Holder 
continuous functions (see Sec. 4.3 and Fig.4.3.5 in [6]). 

Let us now turn to the proof of the theorem. The main ideas and the strat- 
egy of the proof are the following: (i) after having written the equation for the 
invariant manifold, we solve it recursively; (ii) the result of the recursion is conve- 
niently expressed in terms of tree diagrams; (iii) these can be easily shown to give 
rise to a convergent expansion, provided that the large factors e _1 produced under 
iterations by the small rate of contraction of phase space (of order 0(e)) is compen- 
sated by suitable cancellations, stemming from the condition that 7 (x,u>o) = 0. 
All these steps are discussed in detail in the following Sections 3-6. Finally, in 
Section 7, we add some comments on the assumptions made and on possible gen- 
eralizations of our result. 

3 The equation for the invariant manifold 

Let W(x) = u>o + U(x), x(Q) = x and for < t < 2tt write w(t) = wo + 1 + u(x, t), 
with u(x,0 + ) = U(x). Substituting into Eq.(2.1), for all < t < 2ir we find: 



u(x, t) = e g(Sx, wq + t + u(x, t),t). 



(3.1) 



Defining 



7o OM) = g(Sx,w + t,t), 7i (x,t) = —g(Sx,w + t,t), 
G(x,t,u) = g(Sx, Wq + 1 + u,t) — jo(x,t) — uji(x,t), 



(3.2) 



Eq.(3.1) can be rewritten for all < t < 2tt as 




(3.3) 



with 



G(x,t,u) = Y,G p (x,t)u p . 

p>2 



(3.4) 
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Integrating from to t G (0, 2tt], the equation for u(x,t) becomes 

u{x, t) = e^'^Uix) + e £ dr e eV{x ^ ( 7o (x, r) + r, u(x, r))) , (3.5) 



where 



T{x,t,r) = / dr , 7i(x,r / ). (3.6) 



By construction G vanishes to second order in it, hence, if <5 is fixed once and for 
all, one has u(x, 2ir) = e r£ it(x,0) + 0(e 2 ), so that 

Kx,2tt)| < e Fe/2 \u(x,0)\, if - < |u(x,0)| < 5, (3.7) 

provided e is small enough. Therefore the region \u\ < 8 is an attracting set for e 
small (compared to S). 

The condition of invariance for W(x) reads u(x, 2ir) = U(Sx), which gives 

U(Sx) - e £T U(x) = e J** dr e er ^' T ) ( 7o (x,r) + G(x, r, u(x, r))) , (3.8) 
and hence 

oo 

= Y,e £r(k ~ 1) F(S~ k x), 

k= \, (3 - 9) 
F(x) =e J * dTe £r ^ 2n ^( 7o (x,T) + G(x,T,u(x,T))). 



In order to solve the two coupled equations Eqs.(3.5)-(3.8), we find convenient to 
consider the two-parameters (/i, e) equations 

u(x,t) =e l * r(x > t ' 0) U(x)+e f' dre^ r( - x ^ T) (j (x,t) + G(x,t,u(x,t))) , 

J° / / (3-10) 

U(Sx) - e^ T U{x) = e J dr ^frr) ( 7o (x, r) + G(x, r, u(x, r))J , 

first showing existence of a solution analytic in e and Holder-continuous in x, at 
fixed (i, and subsequently proving that fx can be taken to have values fj, > e a little 
beyond e. 

4 Recursion equations 

Eq.(3.10) can be solved recursively as follows. Let 
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u(x, t) = e^^Uix) + £(x, t) (4.1) 

and write £(x, t) = e£i(x, t) + e 2 $,2(x, t) + • • • and U (x) = eU\(x) + e 2 C/ 2 (x) + • • •■ 
Plugging the expansions into Eq.(3.10) and setting A = e Mr , we find 

h(x,t) = f 'dTe^ x ^ 7 o(x,r), h(x,t) = , 
Jo 

r2n k (4-2) 

U x {x) = Y \ k - 1 / dre^ r(5_ x ' 2n ' T ho(S- k x,T), U 2 (x) = 0, 

and, for n > 3, 

ft 00 p 
e„(x,t)=/ drE E e^^G p (x,r) n(e^' T ' 0) ^(x) + e ni (x,t)), 

J ° p=2 ni,...,n p >l i=l 

+ .. .+np- n— 1 

C/ n (Sx) - e» r U n (x) = drJ2 E e" r ^ 27r ^G p (*, r) • (4.3) 



p=2 "I, 

Til + . - . +np=n— 1 

P 



1=1 



which give a recursive definition of £ n ,U n . Note that the second of Eq.(4.3) can 
be solved in a way analogous to Eq.(3.9), so that 



Un (x) = g A"" 1 / dr E E e^ fe *' 2 ^G p (S- fc x, r) • 

fc=l p=2 ni,...,n p >l 

n 1 + ...+n p =n-l (4-4) 

• f[ (e^ 5 - fc ^°)c/ ni (rt,r) + ^(5- fc x,r)). 
i=i 

For the incoming discussion, it is useful to note that by the analyticity of g, there 
exists a constant Co such that 

\G p (x,t)\<C%, e^ x ^<C , (4.5) 

uniformly in x, t, r. Moreover, if we consider the Fourier expansion 

G p (x,r) =E e il/ - x G P (v,r), e^^ = E e^f (i/,t,r), (4 g) 

there exists k > 0, such that 

\G p (u,t)\ < Coe- K M, |f (i/,t, r)| < Coe^l, (4.7) 

uniformly in t, r. 
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5 Tree expansion 



The recursion for £ n , U n can be conveniently represented graphically in terms of 
tree graphs. We start with some basic definitions. A connected graph Q is a 
collection of points (nodes) and lines connecting all of them. We denote with 
V{Q) and L(Q) the set of nodes and the set of lines, respectively. A path between 
two nodes is the minimal subset of L(Q) connecting the two nodes. A graph is 
planar if it can be drawn in a plane without graph lines crossing. 

A tree is a planar graph Q containing no closed loops. One can consider a tree 
Q with a single special node vq: this introduces a natural partial ordering on the 
set of lines and nodes, and one can imagine that each line carries an arrow pointing 
toward the node vo thus establishing a partial order y on the tree. We can add 
an extra (oriented) line £q exiting the special node vo; the added line £q will be 
called the root line and the point r it enters (which is not a node) will be called 
the root of the tree. In this way we obtain a rooted tree 6 defined by V(9) = V(Q) 
and L{6) = L(g)u£ . 

Hence given two nodes v , w € V(6), it is v -< w (or w >- v) if w is on the path 
connecting v to the root line. 

A labeled tree is a rooted tree 6 together with a label function defined on the 
sets L{6) and V{9). 

We shall call equivalent two rooted trees which can be transformed into each 
other by continuously deforming the lines in the plane in such a way that the 
latter do not cross each other (i.e. without destroying the graph structure). We 
can extend the notion of equivalence also to labeled trees, simply by considering 
equivalent two labeled trees if they can be transformed into each other in such a 
way that also the labels match. 

We can identify a line with the nodes it connects; given a line i = wv we say 
that t enters w and exits v , and we shall write also i = l v . Given two comparable 
nodes v and w, with v ~< w, we denote with V(v,w) the path of lines connecting 
w to v , with v and w being included. For any node v -< v o we denote by v' >- v 
the node immediately following v, hence £ v = v'v, and we set v' Q = r. 

Associate with the nodes and lines of any tree 6 some labels, according to the 
following rules. 

With each node v we associate the labels (rj v , k v ), where rj v = 0, 1, while k v = 
if rj v = and k v G N (i.e. k v > 1) if r] v = 1; we call r] v the type label of the node 
v . Furthermore we define 



v£P(v,vo) 

For each node there are p v > entering lines. If p v = we say that v is an 
end-node. If p v > then p v > 2. For p v > 2 let s v be the number of nodes w with 




'V 



(5.1) 
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w' = v such that r] w = 1. 

Represent graphically the node v as a big white ball (white node) if r] v = 1 
and as a small black ball (black node) if r] v = 0. Thus s v is the number of white 
nodes immediately preceding v. 

With each node v we associate a time variable r v € (0, 2tt], with t v i = t, a 
node factor F v = F v (x,t v /,t v ) defined as 

f e /,r(s-^) w , r „) G ^ ((S -/c(-) X)T?) ) e ^r(5- fe ^)x,r„,o) ) % = 0j 

\A^- 1 e^( s "* ( " ) ^> T »)G^( 1 S , - fc ( w )x,r t ,)e^» r ( s "* ( " )ie ' T ''' ) ) % = 1, 

and an integral I v (rj v ), with 

I v (0) = T' dr v , I v (l) = [ * dr v . (5.3) 

JO JO 

With the definitions above, we denote by T n ,u,r) the set of trees with n nodes 
and r) VQ = r]. Then one can check (for instance by induction) that 

U n (x) = Val(z;0), £„(M) = £ Val(x;0), 
6»eT„,i eer„, 

Val(x;0) = JJ I v (i] v )F v (x,t v ,,t v ), (5.4) 

where the integrals are understood to be performed by following the tree ordering 
(i.e. by starting from the end-nodes and moving toward the root). Note that in 
the formula for £ n (x, t) in Eq.(5.4), the variable t v i = t plays no role, as I Vo (l) an d 
F Vo (t,T VQ ) in fact do not depend on t. 

For instance i) and U\{v) are represented graphically as in Figure 1. 



£i(s,t) = - U 1 (x)= - 

Figure 1. Graphical representation of £i(x,t) and ?7i(x). 

Since £2, U2 = the next simplest example is the representation of £3,1/3. For 
instance, say, Us(x) is described in Figure 2. 



u 3 (x) = cf + cf + — a + 






Figure 2. Graphical representation of U^{x). 

In a similar way we can represent graphically the other C/ n 's or £ n 's. For 
example in Figure 3 we represent some contributions to U${x): all the other con- 
tributions are obtained by replacing some white nodes with black nodes and vice 
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versa, and possibly permuting the lines entering the nodes (the only node which 
must be necessarily white is the node vq). 




Figure 3. Graphical representation of some contributions to U§(x). 

In drawing the graphs in Figures 1-3 we have used that p v > 2 whenever p v ^ 
and that 6,2(x,t) and U2(x) vanish identically, as it follows from (4.2). 



6 Convergence 

Given a tree i9 £ 7~ n ,x,r]i f] £ {0, 1}, let be the number of its end-nodes, and 

let -Eo(^) and E\($) be the number of end- nodes of type and 1, respectively. 
Similarly, let N % {&) be the number of "internal" vertices of i? (i.e. the number of 
vertices of # with p v > 0) and iV*(#) be the number of internal vertices of type r\, 
with ri = 0,1. Note that, since p v > 2 whenever p v ^ 0, we have E(i)) > (n + l)/2. 
A key remark is that, thanks to the condition that Jo(x, wq) = 0, the contribution 
associated to every end-node v of type 1 is of order fx, 



[ ^dr^A^-V^ 5 k(v)x > 2n ' T ^j {S- k( - v) x,T v ) < (const.)/i. (6.1) 
Jo 

Using this cancellation and the bounds Eq.(4.5), we find that, for all n € N and 

|ValOz;0)| < ClV £lW II ^""^ ( 6 - 2 ) 



all v e Z 2 , 



v£V(6) 
Vv = l 



for a suitable constant C±, uniformly in x. 

Therefore, in order to obtain a bound on U n (x) we have to estimate the number 
of labeled trees in T n ,i- For each node v with rj v = 1 the corresponding sum over 
k v produces a factor C2/tt _1 . Finally the number of unlabelled trees with n nodes 
is bounded by C3 , with C3 = 4. 

By collecting together the bounds above we find 

\U n (x)\ < C n n~ N ^\ C = C X C 2 C Z . (6.3) 
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Using JV}(0) < AP(tf) =n- E(0) <n-(n + l)/2, we find 

\U n {x)\ <C n ^ n -V' 2 , (6.4) 

which yields a radius of convergence in e given by Eq = C -1 ^ 1 / 2 , hence we can 
take ii = e. 

To study the regularity in x we need to get a similar upper bound for 

\U n {x)-U n {x')\ 

\x-*\P ' /3> °- (6 - 5) 

We proceed as follows. Given a tree •& we introduce an arbitrary ordering of the 
nodes by setting V(i?) = {v±, . . . ,v\yr$\\} and define VL - (i?) = {vi, . . . , v p } and 
V p + = {vp+i, . . . , U|y(^)|}, with the convention that ^o _ (i?) = V^y^(-&) = and 
Vj~ w| (i?) = U +(tf) = V(0). Then we define 

Val±(x;tf) = JJ I v (r] v )F v (x,T v/ ,T v ), (6.6) 

with Valg (x;$) = Valjy^|(x;i?) = 1, and write 

?7 n (x) - = ^ (Val(x;tf) - Val(x';tf)) 

t»eT„,i 

= E E Val;_ 1 (x; 1 ?)Val+. 1 (x';^)- ( 6 - 7 ) 
t)eT„,i p=i 

Each summand in (6.7) can be studied as in the case j3 = 0, with the difference 
that the node factor is F v (x, t v >, t v ) for v € V^^fi), F v (x' , t v i , t v ) for v G ^>+i($) 
and F„(x, r„/, r„) — F v (x' ,t v i ,t v ) for i> = u p . Using Eqs.(5.2) and (4.6), we find 
that 

F„(x,r„/,r„) = e tu ' x F v (u,T v ,T v >) , 

z,eZ 2 (6.8) 
\F v (v,t v ,t v ,)\ < \^- 1 )^C 1 CF +2 e-® sk{v)v \ , 

for suitable constants Co, Ci, independent of r„, r„/ (in deriving the bound on the 
Fourier coefficients we used that F v (x,t v > ,t v ) depends on x through S~ k ^x, see 
Eq.(5.2)). Therefore, 

I F v [x , T v ' , T v ) -P^ (x , 7V , Tj; ) I 

\x-x'\P 

< x (k v -i)8 Vv , l0i0 s v+ 2 j2 \ufe-^\\-\-^ , 
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where A_ is the eigenvalue of the hyperbolic automorphism S with smallest abso- 
lute value; note that, if A+ is the other eigenvalue of S, we have |A_| = |A+| _1 , with 
|A+| > 1. The factor |A_|-^W in Eq.(6.9) is bounded by UveV^y. Vv =i \*+\ Pkv , so 
that 

-Un{x')\ ^ ^ n cy^W FJ A^IA+I^, (6.10) 
' x x I 0eT n ,i «6V(*) 

Vv = l 

for a suitable constant C. We must require that (A+I^A < 1, and in order for the 
sum ^k^id^+f ^) kv to be still proportional to /i -1 , we can take 

P = cfi, < c < |r|/log|A+|, (6.11) 

which gives (3 = ce for fi = e. This concludes the proof of the theorem. 



7 Concluding remarks 

Let us add a few comments on the role of the specific hypotheses in the theorem. 
The assumption that there is wq such that J 2n dt g(x, t, wq + t) is zero has been 
heavily used as well as the assumption that J n dt d w g(x, t, wq + 1) = T, for T < 0. 
Both assumptions can be relaxed into 

(a) J 27r dt g(x, wq + t, t) = eg{x) , for some g(x) , 

(b) J 27r dt d w g{x, w + t,t) < T, for T < 0, 

and the statement of the theorem remains valid (and its proof is essentially un- 
changed) . 

Assumption (a) is strong and it does not allow for perturbations of the form 
eg because the e-dependence has to involve higher powers of e to satisfy (a) unless 
g(x) = (as in the theorem proved here). 

It is unclear whether the results discussed in this paper would hold under the 
weaker assumptions that 

(a') J 2lT dtg(x, u>o + t,t) = g(x), with g{x) with average, 

(b') fo^ dt d w g(x, wo + t,t) = gi(x), with gi{x) with negative average. 

Preliminary computations and numerical evidence suggest that the existence 
of an attractor A = {(x,0,W(x))} for may still be true, possibly "just" for 
a.e.-x and with a non-smooth function W(x): in other words, a similar phase- 
locking phenomenon as the one proved in this paper may still take place, with A 
a fractal rather than a smooth surface. 
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It would be interesting to check that the central Lyapunov exponent in the 
theorem is also expressed by an expansion in tree graphs: the procedure should 
be the same illustrated in the theory of the tangent map for perturbations of the 
simple Anosov map (e.g., A2 in Section 1), see proposition 10.3.1 in [6]. 

Finally a natural conjecture is that the assumption / = in Eq.(2.1) can be 
eliminated as long as / is analytic. However, we do not have specific indications 
of what is the general structure of the attractor in this case. 
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